Signature of elasticity in the Faraday instability 
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We investigate the onset of the Faraday instability in a vertically vibrated wormlike micelle 
solution. In this strongly viscoelastic fluid, the critical acceleration and wavenumber are shown to 
present oscillations as a function of driving frequency and fluid height. This effect, unseen neither 
in simple fluids nor in previous experiments on polymeric fluids, is interpreted in terms of standing 
elastic waves between the disturbed surface and the container bottom. It is shown that the model 
of S. Kumar [Phys. Rev. E, 65, 026305 (2002)] for a viscoelastic fluid accounts qualitatively for our 
experimental observations. Explanations for quantitative discrepancies are proposed, such as the 
influence of the nonlinear rheological behaviour of this complex fluid. 
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Since Faraday's founding work 0, the parametric in- 
stability of a vertically vibrated fluid layer has emerged 
as one of the best candidates to study pattern formation 
and nonlinear dynamics. Above a critical acceleration a c , 
an initially flat and quiescent fluid layer driven vertically 
at frequency / goes unstable and gives way to a pattern 
of surface waves that oscillates at half the driving fre- 
quency with a characteristic wavenumber k c 0,Hlg. As 
the driving acceleration is raised above onset, a series of 
secondary instabilities take place, leading to defect dy- 
namics and eventually to spatio-temporal chaos Q . 

So far, the Faraday instability has been mostly studied 
in simple viscous fluids 0, SHI EE El • Recently, interest 
has grown in the effect of vertical vibrations on a layer 
of complex fluid, both experimentally 0i@ and theoreti- 
cally 0, El ■ Due to their microstructure, complex fluids 
display viscoelastic properties, which may affect classi- 
cal hydrodynamic instabilities [llj . However, previous 
experiments on semidilute polymeric solutions did not 
show significant modifications of the Faraday instability 
0, II] ■ The debate has mainly focused on the existence of 
a harmonic response (at /) instead of the classical sub- 
harmonic response (at f 12), a typical viscoelastic effect 
predicted numerically |fj, |lOj and observed experimen- 
tally at rather low frequencies (/ < 40 Hz) together with 
new types of patterns that compete with each other 

In this Letter, we report onset measurements of the 
Faraday instability in a wormlikc micelle solution, a 
strongly viscoelastic fluid well characterized by a sin- 
gle relaxation time r. The critical acceleration and 
wavenumber are shown to present oscillations as a func- 
tion of driving frequency (see Fig.^). This effect is inter- 
preted in terms of standing elastic waves between the dis- 
turbed surface and the container bottom. It is also shown 
that the finite-depth model of Ref. accounts qualita- 
tively for these experimental observations. Finally, possi- 



ble explanations for discrepancies are discussed and the 
implications of our results in the field of complex fluid 
hydrodynamics are emphasized. 
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FIG. 1: (a) Critical acceleration a c and (b) critical wavenum- 
ber k c vs driving frequency in a 4 % wt. CPCl/NaSal worm- 
like micelle solution (•). Also shown are the calculations 
for the corresponding Maxwell fluid (dotted lines) and when 
adding a Zimm-like term (solid lines). Insets: a c and fc c for 
a Newtonian fluid (78 % glycerol-22 % water mixture) and 
corresponding predictions using the model of Ref. @ with 
density p = 1.19 gem -3 , viscosity r\ = 0.05 Pas, and surface 
tension a = 0.06 Nm" 1 (solid lines). 

Wormlike micelles (also called "living polymers") are 
long, c ylin drical aggregates of surfactant molecules in so- 
lution [111 ] . They spontaneously form under given condi- 
tions of temperature and concentration for a wide range 
of surfactants. Unlike conventional polymers, their size 
is not fixed and they constantly break and recombinc 
under thermal agitation. Their dynamics thus result 
from both a classical reptation motion [J^ and break- 



2 



ing/recombination processes [l3j. This peculiar feature 
leads to an almost perfect Maxwcllian behaviour in the 
small-deformation regime, for which the complex viscos- 
ity reads: t]m{oj) = Gqt/ (l + iuir), where Go is the shear 
modulus, t the relaxation time, and to the pulsation. 

Our working fluid is a wormlike micelle system made 
of cetylpyridinium chloride (CPC1, from Aldrich) and 
sodium salicylate (NaSal, from Acros Organics) dissolved 
in brine (0.5 M NaCl) j3 0- In this study, we fo- 
cus on a 4 % wt. sample with a concentration ratio 
[NaSal]/[CPCl]=0.5 as described in Ref. Linear 
rheological measurements were performed in the cone- 
and-plate geometry using a stress-controlled rheometcr 
(AR 2000N, TA Instruments) with small stress oscilla- 
tions of amplitude 2 Pa (the linear regime extends up to 
20 Pa). As shown in Fig. El the dynamic moduli are well 
accounted for by a Maxwell model G'(to) = — tu^s(r]M(to)) 
and G"{to) = tor) s + w3?(?7a/(w)), with Go = 16 Pa 
and r = 0.44 s, at least for frequencies below 2 Hz. 
i] s = 10~ 3 Pas represents the solvent (brine) contribution 
to the viscosity. 
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FIG. 2: Storage modulus G (to) (o) and loss modulus G"(uj) 
(•) as a function of oscillation frequency u/2ir for a stress 
amplitude of 2 Pa. The dotted line is the best Maxwell fit 
(Go = 16 Pa and r = 0.44 s) with a solvent viscosity rj a — 
10 -3 Pas. The solid line corresponds to the same Maxwell 
model with an added Zimm-like term (see text). 

Our experimental setup is a classical one for studying 
the Faraday instability. A cylindrical container of inner 
diameter 60 mm is filled to a height h = 10 mm under the 
brimful boundary condition [f| and vertically vibrated 
using an electromagnetic shaker (Ling Dynamic Systems 
V406). The fluid temperature is regulated to 21±0.05°C 
by water circulation beneath the container. In order to 
prevent evaporation and surface contamination, the con- 
tainer is sealed by a Plexiglas cover that supports a small 
piezoelectric accelerometer (Endevco 2224C). The signal 
from the accelerometer is fed to a lock-in amplifier (Stan- 
ford Research Systems SR810) that provides the sinu- 
soidal excitation at frequency / to the shaker. The ex- 
periment is illuminated from above by a circle of diodes 
strobed at / or //2. The bottom of the container is made 
of aluminum to provide mechanical rigidity, and anodized 



to prevent light reflection. 

For a given driving frequency, the instability threshold 
is determined by (i) noting the acceleration for which the 
instability first appears, (ii) fully destabilizing the whole 
surface by increasing the acceleration about 5 % above 
onset, and (iii) decreasing the acceleration and noting 
the acceleration for which the instability completely dis- 
appears. In all cases, the two values of the critical accel- 
eration differ by less than 2 %: no significant hysteresis 
is observed. The value of a c shown in Fig. da) is the 
mean of these two accelerations, their difference remain- 
ing always smaller than the marker size. Moreover, the 
surface response was found to be subharmonic (at //2) 
over the whole range of investigated frequencies / = 25- 
90 Hz. For each frequency, the critical wavenumber k c is 
estimated from a picture of the fully destabilized surface. 

Figure H constitutes our main result. The critical ac- 
celeration and wavenumber are not monotonic but rather 
oscillate with the frequency, a maximum in a c corre- 
sponding to a large drop of k c . To our knowledge, such 
a marked effect has not been reported in previous ex- 
periments on complex fluids. In Newtonian fluids, non- 
monotonic a c or k c are associated to lateral boundary ef- 
fects (as only an integer or semi-integer number of wave- 
lengths may fit in the cell when the brimful condition 
is used) 0, [Rj. In order to rule out such an inter- 
pretation, we repeated the experiment with a Newto- 
nian fluid of critical acceleration and wavenumber sim- 
ilar to our wormlike micelle solution. The insets in Fig.^ 
show monotonic behaviours, in perfect agreement with 
the model of Ref. [|| for a laterally unbounded viscous 
fluid, as already found in other Newtonian fluids This 
allows us to consider our results in terms of a laterally un- 
bounded fluid and confirms the viscoelastic nature of the 
oscillations seen in a c and k c for the wormlike micelles. 

More precisely, we propose to interpret these oscilla- 
tions as a coupling between the disturbed surface and 
clastic waves reflected at the container bottom. For a 
Maxwell fluid in the lot 3> 1 limit, which is always ver- 
ified in our experiments, the velocity c and attenuation 
coefficient a of shear waves are given by c = \jG§j p and 
a = 1/(2ct), where p is the fluid density. Above on- 
set, surface disturbances generate elastic waves that may 
form a standing wave in the container height h. When 
2h is a semi-integer multiple of the wavelength 2c// [l7j . 
i.e. / = (n + 1/2) /i with /i = Go/p/h and n an inte- 
ger, constructive interferences form at the surface. In this 
case, elastic waves are amplified and surface disturbances 
are promoted so that the critical acceleration is expected 
to be lower than in the absence of elastic effects. On 
the other hand, when / = nfx, destructive interferences 
should lead to an increase of a c , hence the oscillations of 
Fig.^a). Moreover, at each maximum of a c , the number 
of elastic modes in the vertical direction increases. Since 
the vertical wavenumber is linked to the horizontal one 
through the incompressibility condition 0,0, jumps sim- 
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ilar to those of Fig. ^b) are expected even in a laterally 
unbounded fluid. 

In order to further check the above interpretation, 
we performed the finite-depth numerical calculations of 
Ref. for a Maxwell fluid with the rheological parame- 
ters inferred from Fig. [2 The two remaining parameters 
used in the calculations, namely the density and the sur- 
face tension of the fluid, were measured independently 
yielding p = 1.0 gcm~ 3 and a = 0.03 Nm" 1 respectively. 
By varying h (Go resp.) for a fixed Go (h resp.), we show 
in Fig. that large oscillations also show up numerically. 
Moreover, they perfectly agree with the simple relation 
fn = n \/ Go/ p/h suggested above. The interpretation in 
terms of clastic waves reflected at the container bottom is 
thus confirmed. This finite-depth viscoelastic effect went 
unseen in previous numerical studies that focused on the 
question of harmonic vs subharmonic response 0, Hoj . 
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Finally, both our simple analysis and the numerical re- 
sults of Fig. EJa) predict that a c and k c should oscillate 
when the fluid height h is varied at a fixed frequency. Fig- 
urc^shows that such oscillations indeed take place in the 
experiment and that the Maxwell model corrected by a 
Zimm-like term provides a good description of a c . Again, 
these findings are strikingly different from the monotonic 
decay of a c vs h in Newtonian fluids H . 
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FIG. 3: Numerical calculations using the model of Ref. 
with t = 0.44 s, p = 1.0 gem" 3 , and a - 0.03 Nm" 1 . (a) a c 
vs h and / at fixed Go = 16 Pa. (b) a c vs Go and / at fixed 
h = 10 mm. White dots are the lines f n = n\l Go/ p/h with 
n = 5 and 7 and coincide with maxima of a c . 

Note, however, that these numerical calculations based 
on the Maxwell model do not allow any quantitative com- 
parison with the experimental a c and k c (see dotted lines 
in Fig. nj . This is most probably because our fluid does 
not remain purely Maxwellian at frequencies higher than 
2 Hz. Indeed, although our rheological measurements are 
limited to 10 Hz, an upturn of G"(oj) is clearly visible on 
the data of Fig. |2 Such an upturn is usually attributed 
to the presence of Rouse- like or Zimm-like motion at high 
frequencies . If we model this behaviour by adding 

a Zimm relaxation term to the previous Maxwell model, 



rj(uj) =r) s + 
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a very good fit of both G'(w) and G"(w) is obtained with 
a z = 0.12 Pas 2/3 (see solid lines in Fig.0). 

Using Eq. JQ) in the numerical calculations leads to 
computed curves (shown as solid lines in Fig. ^| that 
are much closer to the experimental measurements, al- 
though the attenuation of the oscillations with increasing 
frequency is now too strong (see discussion below). 



FIG. 4: (a) Critical acceleration a c and (b) critical wavenum- 
ber k c vs fluid height measured at 50 Hz (□) and at 70 Hz 
(•). The solid lines correspond to the Maxwell model with an 
added Zimm-like term (see text). 

Let us now discuss the present experimental and nu- 
merical results. Previous experiments on polymeric solu- 
tions did not show any oscillations in the onset measure- 
ments 0, HI • This can be explained by considering the 
attenuation coefficient of the elastic waves a. If a > 1/h, 
no standing wave, and thus no oscillation, is observed 
due to a too strong attenuation. In our wormlikc micelle 
solution, the storage modulus G'{ui) remains always at 
least three times larger than the loss modulus G"(w) and 
ah ~ 0.09, so that oscillations are observed, whereas pre- 
vious experiments were performed either in the capillary 
regime or in the overdamped regime || where elastic 
waves may not be excited [ijj . 

Moreover, even though the oscillations in a c are rather 
well reproduced by the Maxwell model with a Zimm cor- 
rection, the agreement with the experimental values of 
k c remains only qualitative (see Figs.QJb) and^Jb)). A 
Rouse correction was also tried with similar results. It 
may be argued that the discrepancy is due to an inad- 
equate rheological model. Indeed, standard rheometers 
are limited to about 15 Hz, so that rheological data had 
to be extrapolated to higher frequencies 

However, we believe that the lack of agreement points 
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to more fundamental questions about the validity of the 
numerical approach. Indeed, besides lateral finite-size 
effects, which we ruled out at the beginning of this study, 
at least two important physical effects are not taken into 
account in the model of Ref. : interfacial rheology and 
nonlinear bulk rheology. 

First, in the numerical calculations, the air-fluid in- 
terface is modelled as purely elastic with surface tension 
(j, whereas in a complex fluid, a can depend on the fre- 
quency. More importantly, the interfacial viscoelasticity 
may significantly affect the boundary conditions |20l |. 

Second, the model of Ref. relies on linear viscoelas- 
ticity. Indeed, since the instability occurs from a qui- 
escent state, strain is expected to be "small." Experi- 
mentally, the smallest detectable surface deformation is 
£ ~ 10 fjxa. The corresponding strain rate 7 may be es- 
timated by 7 ~ (uk ~ 10 s _1 at / = 60 Hz p. Since 
the zero-shear viscosity of our fluid is 770 = Gqt ~ 7 Pas, 
7 ~ 10 s -1 corresponds to a shear stress of about 70 Pa, 
which is far into the nonlinear regime. Thus, the only 
knowledge of the linear viscoelastic moduli G' and G" 
may not be sufficient to fully account for the experiments 
and nonlinear effects could provide an explanation for the 
observed discrepancies. 

More precisely, wormlike micelles are known to align 
under shear, leading to a shear-induced ncmatic state 
This strongly shear-thinning isotropic-to-nematic 
transition occurs for 7 ~ 0.5-50 s _1 depending on tem- 
perature and concentration |l5j . In a Faraday wave 
pattern, shear is localized between the crests and the 
troughs. Thus, above onset, aligned micelles are expected 
to coexist with the isotropic state and to spatially mimic 
the surface pattern. Finally, since the destabilized state 
is far from a pure shear flow, one may also wonder about 
the influence of nonlinear extensional rheology 

To conclude, the present study reveals a strong signa- 
ture of elasticity on the onset of the Faraday instability 
in a wormlike micelle solution. This original effect re- 
sults from standing elastic waves in the container height 
and should be very general in highly viscoelastic fluids. 
It shows that, as suggested earlier the character- 

istic relaxation time of a complex fluid may couple to 
the excitation period and/or to memory effects, leading 
to new temporal and spatial behaviours under vertical 
vibrations. Linear viscoelasticity allows for a good qual- 
itative agreement with the model of Ref. but we be- 
lieve that surface rheology and/or nonlinear effects are 
also significant. This raises new challenges for theory 
and modelling of hydrodynamic instabilities in complex 
fluids. In particular, we point out that the presence of a 
shear-induced isotropic-to-nematic transition could play 
a major role in the pattern selection. This last remark 
has to be related to some very recent experimental find- 



ings on vertically vibrated shear-thickening suspensions, 
where a rich variety of patterns such as holes and fingers 
was observed after a finite perturbation was applied to 
the surface |2 ll |. Such an unusual behaviour was linked 
to the nonlinear rheological properties of the fluid. Our 
results on wormlike micelles provide another example of 
a striking effect induced by the microstructure of a com- 
plex fluid on a classical instability. 
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